Propagation of solitons of the magnetization in magnetic nano-particle arrays 
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It is clarified for tlie first time tliat solitons originating from 
the dipolar interaction in ferromagnetic nano-particle arrays 
are stably created. The characteristics can be well controlled 
by the strength of the dipolar interaction between particles 
and the shape anisotropy of the particle. The soliton can 
propagate from a particle to a neighbor particle at a clock 
frequency even faster than 100 GHz using materials with a 
large magnetization. Such arrays of nano-particles might be 
feasible in an application as a signal transmission line. 



be feasible in an application as a signal transmission line. 
We examine a one-dimensional array of particles in the 
x-direction. Each particle has a magnetic moment Mi 
{\Mi\ ^ M). The dipolar interaction and anisotropy en- 
ergy is modeled as 
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The recently developed nanofabrication techniques 
make it possible to fabricate ferromagnetic particles to 
a length of 20-30 nanometers.^ Since the size of such 
nano-particles is small enough to be a single domain, 
the magnetization is homogeneous over a particle and 
can be described by a magnetic moment. These parti- 
cles interact with each other due to a dipolar interac- 
tion of the magnetic moments. In such a nano-particle, 
the magnetic anisotropy energy can be dominated by 
the particle's shape rather than the magnetocrystalline 
anisotropy. Therefore, the magnetic anisotropy energy 
can be well controlled by changing the shape of the parti- 
cle. When the distance between particles is short enough 
that dipolar interaction becomes comparable with the 
anisotropy energy, the direction of the magnetization of 
each particle is determined in order to minimize both en- 
ergies, and the direction of the magnetization moves in a 
corrective manner among the particles. 

It is known that the magnetic domain wall in bulk mag^. 
netic materials is a kind of solitons of nonlinear waves.cTcl 
Usually, the soliton originates from the exchange inter- 
action between spins. This letter clarified for the first 
time that, in a one-dimensional array of nano-particles, 
solitons originating from the dipolar interaction between 
particles are present, and that the characteristics can be 
well controlled by the strength of the dipolar interaction 
and the shape anisotropy. The great advantage of con- 
trol is quite a contrast to the case of the soliton of the 
domain walls in bulk materials. 

In this letter, we analyze the characteristics of the soli- 
tons in arrays of ferromagnetic nano-particles by vary- 
ing the experimentally controllable parameter: the shape 
anisotropy energy. When the damping is small, the soli- 
ton propagates from a particle to a neighbor particle at 
50 GHz, even faster than 100 GHz using materials with a 
large magnetization. Such arrays of nano-particles might 



where Ky and K^ describe the shape anisotropy. The 
dipolar interaction is characterized by J — hqM'^ / {^tto?) . 
In this letter, we consider the case of K^ ^ J and 
Kz ^ Ky, and thus the x — y plane is an easy plane 
for the magnetization. Further, we consider the case 
where the thermal fluctuation in the direction of Mi, 
which frequently leads to the super-paramagnetism,l3 is 
well-suppressed by Ky and/or J . As shown below, the 
condition can be sufficiently satisfied since J can exceed 
10, ODD K even for particles of a size of 20 nm, when these 
are aligned closely to each other. 

The magnetisation of each nano-particle obeys the 
Bloch equationu 
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where u = g^isl^ is a gyromagnetic constant, and a 
describes the dumping due to energy dissipation, such 
as magnetoelastic dissipation and eddy current loss in 
metallic ferromagnets. In nano-particles, however, it has 
been shown that the amount gf these dissipation mech- 
anisms are negligibly smallfl'EI In fact, considering the 
eddy current loss, a ~ 10"'' for Fe-particles of a size of 
20 nm, for which the soliton can propagate beyond 1,000 
sites. Therefore, we exclusively consider the case of a = 
in this letter. 

Before discussing the results of the numerical simula- 
tions of Eq. (0), it is worth considering the continuum 
limit. There are two types of configurations of Ali giv- 
ing the minimum total energy depending on Ky. When 
only the nearest-neighbor interaction is considered for 
simplicity, the configuration is type I [Fig. |l] (a)] and 
type II [Fig. |^ (b)] for Ky < J and Ky > J, respectively. 
In both cases, the Bloch equatioa (||) is approximated to 
the sine-Gordon (SG) equationiQ 
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FIG. 1. The configuration of the magnetization giving a 
minimum total energy, (a) type I for Ky < J and (b) type 
II for Ky > J. For both cases, an example of the soliton is 
shown below each figure. 



where c^ = 3Ja^K,y/AP, X^ = 3Ja^/{A\Ky - J\), and 
0'{x, t)/2 — 9{x, t) is the azimuthal angle of M{x) in the 
continuum limit of Mi. Note that the definitions of Oi 
for even- and odd-sites are different as shown in Fig. |l|. 
Consequently, the solution in the form of the soliton is 
present for both cases of type I and II, whose examples 
are shown below Figs. H (a) and (b). The solution has 
a velocity (w) as a parameter, and its wavelength is / = 
A-\/l — v^ jt? . The maximum velocity is limited to c. 

When the long range part of the interaction is taken 
into account, \Ky—J\ in A^ is replaced by liCy— 5(^(3) J/4| k, 
\Ky — \.hJ\ with C(x) being Riemann's zeta function. To 
gain a large value of J, however, it is better to align 
particles closely in such a way that the distance be- 
tween particles (a) is comparable with the diameter of 
the dot (d), where the nearest-neighbor interaction be- 
comes dominant. Therefore, we only take into account 
the nearest-neighbor interaction in Eq. (|l|) hereafter, al- 
though the actual form of the interaction may slightly 
deviate from that in Eq. (||) for a^d. 

Figure shows the time-evolution of the center coor- 
dinate of the soliton obtained from the numerical calcu- 
lations of Eq. (g) for various values of Ky. As an initial 
condition for i = 0, we choose a solution of SG-equation 
(^ with an initial velocity v = 0.2c. For both cases of 
type I and II, the soliton stably propagates keeping the 
initial velocity when Ky is close to J. The velocity begins 
to decrease when jiiTj^— J| becomes large, and the amount 
of the decrease is more significant for the case of type II 
than type I. The soliton for type II stops after the prop- 
agation beyond 30 sites for Ky = 1.2 J, while the soliton 
for type I can propagate even for Ky = 0. Therefore, the 
configuration of type I with Ky close to J is suitable for 
the stable propagation of the soliton. 

We performed the same calculations for the type I case 
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FIG. 2. The time-evolution of the center coordinate of the 
soliton for various values of Ky . Solid and dotted lines are for 
type I and type II, respectively. A gray zone shows the range 
(200 <x< 400) where the average velocity plotted in Fig. Q 
(a) is obtained. 



for various values of the initial velocity, and obtained 
the average velocity (v) in the range, for example, of 
200 < a; < 400. The results are shown in Fig. | (a). The 
average velocity v keeps the initial velocity for smaller ve- 
locities. This implies that the solution of the continuum 
limit is appropriate for this range of the velocity. With 
the increase of the initial velocity, however, the effects 
of the discreteness of the array become important, and 
V begins to saturate and comes to depend little on the 
initial velocity. With the further increase of the initial 
velocity, v for Ky >0.3J increases suddenly with a finite 
jump, and v takes an almost completely fixed maximum 
value. It is interesting to note that the corresponding 
wavelength of the soliton is always close to the distance 
between particles (a) independent of Ky, while the corre- 
sponding fixed value of v depends on Ky. In the case of 
Ky < O.U, although the jump becomes very broad and 
the velocity is not fixed so rigorously, similar behavior 
also remains. It should be noted that the sudden jump 
and fixing in the velocity is absent in the case of type II. 
According to the numerical results, the maximum of 
the average velocity (vc) and the corresponding wave- 
length (Ic) satisfies an approximate relation similar to 
the continuum limit: 
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as shown in the inset of Fig. ||. Therefore, due to the the 
effects of the discreteness of the array, the wavelength is 
shifted by a constant, and c and A are renormalized si- 
multaneously, although the amount of the shift and the 
renormalization may depend on the distance of the soli- 
ton propagation. The corresponding wavelength Ic is al- 
ways close to a as discussed above. Therefore, the max- 
imum velocity of the soliton is approximately given by 
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FIG. 3. (a) The average velocity obtained for the range 
200 < X < 400 as a function of initial velocities. The inset 
shows the relation between the maximum of the average ve- 
locity (vc) and the corresponding wavelength {Ic). (b) The 
average velocity of the soliton generated by flipping Mq as a 
function of the flipping time. The average is performed in the 
range of 200<a;<400. 



an Fe-particle of the above size, io/\/3ifz/>/~2 x 10 ^^ 
sec (^ 50 GHz) assuming K^j J = 5. The maximum 
velocity in the continuum limit is c = ^JiK^/J (a / to) ^ 
1 X 10"^ m/s for a = 20 nm. When the particles are 
fabricated to have a shape anisotropy of Ky ~0.7J, the 
actual velocity is about 70% of c according to Eq. (||). 
Since the clock frequency shown above is proportional to 
the magnetization, it might be able to exceed 100 GHz 
for materials with a large magnetization. 

To conclude, we have shown for the first time that 
a soliton excitation mode is present in magnetic nano- 
particle arrays with the dipolar interactions for both con- 
figurations of the magnetization (type I and type H) giv- 
ing the minimum total energy depending on the shape 
anisotropy. The soliton mode in the type I configuration 
exists in a wider range of the shape anisotropy than that 
in the type II configuration. The solitons with the max- 
imum velocity, which depend on the shape anisotropy, 
always have a wavelength close to the distance between 
particles. These solitons are stable independent of the 
initial conditions, and they can be generated by fast fiip- 
ping of the magnetization at the edge of the array. The 
soliton can propagate from a particle to a neighbor par- 
ticle at a clock frequency even faster than 100 GHz, and 
may be feasible in an application as a signal transmission 
line. 

Finally, it should be noted that the dominant mech- 
anism of the dumping in the nano-particle arrays and 
the effects of the disorder in the alignment of the parti- 
cles and its shape are not clarified yet. Further intensive 
studies will be necessary on these effects. 

We would like to thank Takeshi Honda, Jun-ichi Fujita, 
and Toshio Baba for helpful discussions. 



It is possible to generate similar solitons when the mag- 
netization at the edge of the array (Mq) is simply rotated 
and flipped by an external force. In fact, the average ve- 
locity of the solitons generated in such a way is plotted in 
Fig. ^ (b) as a function of the flipping time. When flip- 
ping is faster than M/Jv, the average velocity is fixed 
without regard to the flipping time. The corresponding 
wavelength of the soliton is also flxed to a, and the veloc- 
ity is approximately given by Eq. (M) again. For slower 
flipping time (>3M/Jj/), the average velocity gradually 
decreases with the increase of the flipping time, and we 
have checked that the soliton can be generated even for 
much slower flipping time {^lOOM/Ji/). 

Finally, we mention the actual values of parameters 
for an experimental setup. It is better to align particles 
closely (a ~ d) in order to gain a large J as discussed 
before. Since the particle is a single domain, M (x (P, 
and thus J oc M'^/a^ oc (P. Therefore, J decreases with 
the decrease of the size of the particle. For a cylinder of 
Fe with d = 20 nm and a height of 10 nm, however, J 
is as large as J '--^ 2 x 10* K. The soliton p ropagat es to 
the nearest- neighbor site in a time of tt^j ^J'iKzj J with 
tQ = M/{Jv)^ which does not depend too much on d. For 
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